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An exact expression is derived for the (w,p) = thermal correlator of shear stress in 
SU(A^c) lattice gauge theory. I remove a logarithmic divergence by taking a suitable linear 
combination of the shear correlator and the correlator of the energy density. The operator 

o ■ 

p\J ' product expansion shows that the same linear combination has a finite limit when cj — > cxd. 

^^' It follows that the vacuum-subtracted shear spectral function vanishes at large frequencies at 

least as fast as aj. (oj) and obeys a sum rule. The trace anomaly makes a potential contribution 



in 



to the spectral sum rule which remains to be fully calculated, but which I estimate to be 
numerically small for T > 3Tc. By contrast with the bulk channel, the shear channel spectral 
density is then overall enhanced as compared to the spectral density in vacuo. 

PACS numbers: 12.38. Gc, 12.38. Mh, 25.75.-q 



£^ 

>► ■ I. INTRODUCTION 

\o 

00 

^*^ I The shear viscosity r/ is a quantity that universally characterizes the relaxation of a fluid towards 

l/^ ■ equilibrium. Indeed rjc? /[kB{e + p)] is the diffusion coefficient of transverse momentum, and this 

o ■ 

(•^ \ process is in many situations occurring in Nature the dominant mechanism of dissipation (e is 

the energy density and p the pressure of the system). The dimensionless quantity r]/{hs) can be 

jv>( I thought of as the ratio of a transport time scale r//[A;^(e + p)] to the thermal time scale h/{kBT) 

C^ ■ and is therefore a good measure of the ability of a fluid to flow [l|, I4] . 

In the past years, the phenomenology of the RHIC heavy ion experiments has provided an upper 
bound on the ratio rj/ s of hot quark matter (see 3] for a review; from here on we set h, c and ks 
to unity). In parallel to this, the thermodynamic properties of QCD in the range of temperatures 
100 < T/MeV < 700 are the subject of ongoing Monte-Carlo simulations on a space-time lattice [J- 
^. In this computational approach, access to the near-equilibrium properties, such as the shear 
viscosity, is limited because lattice QCD employs the Euclidean formulation of thermal field theory. 
Real-time properties can thus only be determined by analytic continuation, a numerically ill-posed 
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problem (see e.g. [7|, [g] section 5). 

Through the Kubo formula, the shear viscosity is related to the low-frequency part of the 
spectral density, rj = vrlinit^^^o ^12,12/"^ ([9] and Refs. therein). Here and in the following we 
denote by TTPf^u,p(7 the imaginary part of the retarded correlator of the energy-momentum tensor 
components T^u and Tp^^. Some constraints on the thermal behavior of the spectral function 
have been obtained by computing the Euclidean correlator of the stress-energy tensor on the 
lattice 



101, 



ll[ |. Because the spectral density is related by an integral transform to the Euclidean 



correlator, one is automatically led to studying the spectral density over the whole semi- axis of 
frequencies. For that reason, it is helpful to have global constraints on the spectral densities, such 
as sum rules [l^, ll3| . Very recently, I determined some of the gross features of the bulk-channel 
spectral density by combining Euclidean correlation functions with a spectral sum rule [ij]. The 
goal of this paper is to take the steps necessary to apply the same strategy to the shear channel. 

Romatschke and Son jlSf] obtained a shear sum rule for conformal field theories, and proposed 
a modified sum rule that takes into account the trace anomaly, T^^. Here I give a derivation of the 
shear sum rule for the SU(A''c) gauge theory in the lattice regularization. In the continuum, the 
shear sum rule reads 

/oo J 2 
[pi2.i2(w,P = 0,T)-/>i2,i2(w,P = 0,0)] = -e(T)- lim AG{uj,T), (1) 

where AG{uj,T) is defined by Eq. (j3ip : I will show by explicit calculation that its a; —t- 00 limit is 
a finite quantity, and that it is proportional to e — 3p. The right-hand side of Eq. ([T]) is therefore 
consistent with the form of the sum rule given by Romatschke and Son. Although the calculation of 
the coefficient remains to be completed, I estimate its practical importance for the determination of 
the shear viscosity in the range of temperatures explored at the LHC and find it to represent at most 
a 5% correction to the energy-density term on the right-hand side of Eq. ([T]). As a consequence, 
the right-hand side is positive at those temperatures, indicating that p[uj)/uj is overall enhanced 
relative to the vacuum spectral density. This is in contrast with the bulk sum rule, which indicates 
an overall depletion of the bulk spectral weight throughout the deconfined phase. 

I start by deriving an expression for the a; = p = Euclidean shear-channel correlator in section 
(jTI|) . Next I convert this identity into a spectral sum rule for the corresponding spectral density 
(Sec. inij) . The contact term that has to be subtracted before using the spectral representation is 
then shown to be finite in section (HYj). I show that part of the finite term can be obtained from 
existing results, while the other remains to be calculated. The numerical importance of the contact 
term is estimated on the basis of this partial result. I make some final comments in section (jV|l . 



II. DERIVATION OF THE LATTICE SUM RULE 

I start by introducing the essential notation, which follows closely [la ]. The reader is referred 
to that paper for unexplained notation. The stress tensor in the SU(A'c) gauge theory reads 

9ix) = m/c^g) F^^{x)F^^{x) 9^,{x) = i v^;.^;. - Ko^Ka- (2) 

The beta-function is defined by qdg/dq = /3(g) = —g^{bo + ...) and 60 = ll-^c/(3(47r)^), 
61 = 34A''^/(3(47r)^) in the SU(A''c) pure gauge theory. If (. . .)t denotes the thermal average 
at temperature T, 

e-W={e)T-{e)o = {e)T-o, e + P=|(0oo)T. (3) 

I will be considering the Wilson discretization of the SU(A''c) gauge theory on an anisotropic lattice, 

Sg = J2f^-S-(^)+^rSr{x), (4) 

X 

where Sa and St are respectively the sum of spatial and temporal plaquettes. The bare gauge 
coupling is given by g^ = 2Nc/y/P7pT- The time direction is discretized more finely by a factor ^, 
called the (renormalized) anisotropy, than the spatial directions. I use the following discretizations, 

r' Q{x) = Z+(/3^, PT)Sa + Z+{/3^,^T)Sr (5) 

r' eoo(x) = z-(/3^, Pt)s^ - z;{p„,Pr)Sr ■ (6) 

The factor Z^^ are such that (V^. Qqo{x))^{J d^x0oo(a;)) in the continuum limit. 

In the earlier publication [}ld|, I derived lattice sum rules for the (a;,p) = two-point functions 
of 6 and ^oo- In these sum rules, the functions of bare lattice parameters A,]q and Aqq appear, which 
were defined as 

Af„(/3„A)4(|--^)(z-±2-,. (7) 

In the appendix of [16|], I gave the values of the latter in the limit g^ -^ 0, corresponding to a lattice 
spacing exponentially smaller than the confinement scale. These values are wrong. The derivation 
of these values was based on the erroneous idea that the operator ^oo does not mix with the unit 
operator even on the anisotropic lattice. I presently correct this mistake. 



A. Perturbative determination of Anr 



00 



The parameters (/3(j,/?t) appearing in the lattice action are related to the renormalized param- 
eters (flo-,^) in a well-defined manner (up to O(a^) ambiguities; the spatial lattice spacing Co- is 
measured in units of a renormalized quantity such as a glueball mass) . It is convenient to exchange 
the 'bare' variables (/3o-,/?t) for {l3,S,o), where 



^0 = V(^r/(3a, 



and 



p = ^/p;^. 



(8) 



The advantage of the new variables is that at the classical level, there is no difference between the 
bare anisotropy ^o a-nd the renormalized anisotropy ^. Since the spatial and temporal plaquettes 
respectively yield a^-B^ and a^a'^E'^, it is not hard to guess the form of Z~^, 



One then easily finds 



z- = ^ 



/? /3'/' 



3/2 ' 



z: 






(9) 
(10) 



'00 



A, 



00 



3 + 0(5o'), 
l + 0(5o')- 



(11) 
(12) 



To determine the quantum corrections, a little more work is needed. First, it was shown in [16l | 
that at the isotropic point, the renormalization factors Z~^ are given by the renormalization of the 
anisotropy. 



Z- ^=' Z- ^=' Zigo) 



di 



(13) 



?=i 



The factors Aqq evaluated at the isotropic point can be written as 



«=i 



'00 



^(.0) + ^ 



9^/3(a^,g) 1 d^/3(fl<.,0 
a(loge)2 16 d(loga)2 



A+ ^= 2 + Z{go) 



Z{go) 






,4 j„-2 



+ 



«=1 



2Z{go)dloga dgl 



(14) 
(15) 



The relation between the bare and renormalized anistropy is known to one-loop order [17]. Indeed, 
two functions c^j^riC) were computed there ([17[], Eq. (2.24)-(2.25)), from which the coefficients 



dc, 



dri€=i ''"^ <- 



cPc 



d( 



f-\^=i can be obtained. The derivatives of the bare with respect to 



the renormalized anisotropy can be expressed in terms of these coefficients, 



Z{go 



de 



l-i(c;-c;)<7o^ + 0(g, 



0)^ 



«=i 



?=i 



-9lK-c'^ + W^-c^)]+Oi9'o)- 



For instance, numericahy [l 



4 - 4 = %J- • 0.146711 - Nc ■ 0.019228 . 



Inspecting Eq. ([H]), we see that the square bracket is at least 0{gQ), and therefore 



Similarly, 



Aoo(5o) = Z{go) + 0{g^o) = 1 " ^(4 - c^)g'o + 0(<7o') • 



Koiao) = 3 + 1 [c'^ - c^^ + c:, - 4] + 0(5o^) . 



These results will be used in sections (IIIII) and (IIVI). 



(16) 
(17) 



(18) 



(19) 



(20) 



B. Sum Rule for 9ii 



In order to be sensitive to correlations of shear stress, we need to derive a sum rule for (say) the 
two-point function of the operator ^n. For that purpose we will need to work on an anisotropic 
lattice, and will take the isotropic limit at the end. To get at the ^n correlator, the anisotropy has 
to be in the 1-direction. 

Since it is conventional to associate the time-direction with the direction where the lattice 
spacing is smaller by a factor ^, I temporarily interchange the labels of the x-direction 1 and the 
time-direction 0. I thus consider a lattice of dimensions ooq ^ i^a ^ 002 x 0x33) and calculate the 
correlator of ^oo in this new coordinate system. At the end I will restore the normal labels, at what 
point the direction with N^ lattice points will play the role of the Matsubara cycle. 

The procedure now closely follows 16|. First, consider a generic renormalization-group invariant 



quantity f^aa,^, T), which is obtained as the continuum limit of a function of the bare parameters, 
F(l3fr, Pr, Nfj). Expressing the independence of f{aa,(,,T) on Oo- and S,, I obtain respectively 



L 



dl 
dL 




d log aa d log Oct 

d log 5 d log 5 




(21) 



At the symmetric point ^^ = 1, the determinant of the matrix is 2/3Z(/3) ^^j^^ . By taking a suitable 
hnear combination, one finds 

This equation holds up to O(a^) effects. 
I now apply Eq. (j22]) to the case of 

f{L) ^ L^iOoo) , (23) 

which is the continuum limit of the lattice expectation value 

F{l3^,Pr,N^) = Ki ((0oo)l - (eoo)(L=oo)) • (24) 

Since we are on the anisotropic lattice, it is necessary to perform the subtraction on the infinite 
lattice (L = c«) in order to remove the mixing with the unit operator. After manipulations entirely 
similar to those performed in [16], I get 

-dK'Wj'^^"" iV^o"o/5^(/5)(^+)^ + <Ao'o(0oo)^ - K E^0oo(^)eoo(O))^ . 

(25) 
Returning to the normal coordinate system, where the compact direction is time and / is given 
by —4(6 + p) at the temperature T = 1/L, the lattice sum rule (j22p can be written as 

«-^E. en(x)Gn(0))T-o - ^^^^\f^^^\ e - 3p) = (1 - iA+o)(e + p) + ^.T^dr^. (26) 

d log a — I 

For completeness I reproduce the sum rule for the 0oo correlator obtained in [l6|, 

a-'{E. eoo(x)eoo(0))T-o - ^^(^o(/?) (^ _ 3^) ^ 3^+ (^)(^ + p) + (!) V^9t^. (27) 

dloga 

It is noteworthy that the difference of two-point functions 
«"'E.(!0ii(^)0ii(O) - fGoo(x)eoo(0))^_^ 

= a-^E.(i(0ii-022)(x)(eii-e22)(o)-|eoo(a:)eoo(o)\ (28) 



= !(l-A+o)(e+p)-iT^9T^ (29) 

is UV-finite. In writing Eq. ()28p I have used the elementary identity 

E^ (eoo(x)eoo(o)) = E. (3 eii(x)eii(o) + 6 611622(0)) . (30) 

A potential application of Eq. (j29p is a determination of the renormalization factor Z^go), since it 
appears quadratically on the left-hand side and only linearly on the right-hand side of the equation. 
However it is probably in the continuum limit that this relation is most useful, as described in the 
next section. 



III. SUM RULE IN THE CONTINUUM AND DISPERSION RELATION 

We are now ready to write a sum rule for the correlator of shear stress, T12. I define 

AG(a;,r) ^ I d^xe^'^^«(i(rii(x) -r22(x))(rn(0) -T22(0)) - |roo(x)roo(0))^_^. (31) 

The lattice sum rule (j29p shows that AG(0,T) is finite. In the next section, I will show that also 
linitj^oo AG(a;, T) is UV-finite. I therefore write the dispersion relation for this correlator as (for 
the general method of deriving sum rules, see [15|, Il8| ) 

/oo ^ 
— (iApn,n(^,r) - iApn,22(w,T) - |Apoo,oo(w,T)) 
-00 "^ 

Api2,i2(w,T)-|Tc^. (32) 



In the last step, I have used rotation symmetry, which implies pi2,i2 = ^Pii,ii ~ ^Pn,22 (the spatial 
momentum p is set to zero throughout this section). I have also used the fact that poo,oo('^)^) = 
TcvUj6{uj). 

On the other hand, AG(0) can be rewritten identically as 

AG(0,r) = j d^x(fen(:E)0n(O) - f ^00(^)^00(0) - ^0(x)0(O) - i0oo(a;)e(O))^__^ . (33) 

In view of Eq. (j29p and in view of the (continuum) sum rules [16|] 



{Jd'xeix)emT-o = T'dT—;^ (34) 



(/ d'x e{x)eoo{0))T = \T^dT^-±l ' (35) 

we have, in infinite volume and in the continuum limit, 

AG(0,r) = ^(e-rc,). (36) 

To reach this expression, we have used the thermodynamic relations Cj, = ^, s = ^, and Ts = 
e + p, as well as Eq. (fTT|) . Now combining Eq. ([32]) and ([36]) . we reach Eq. ([1]) by recalling that, in 
infinite volume, 

/ d^x(^l{Tn{x)-T22{x)){Tu{0)-T22m) = j d? x {Ti2{x) Ti2m . (37) 

In the free theory, we can check the coefficient of the e = ^^ T^ term in Eq. ([T]) {d^ = N"^ — !)■ 
Using 



)(47r)2tanhw/4r ^isJ 
the coefficient of the energy density on the right-hand side of the sum rule is correctly reproduced. 



^^•^^ = m(LrSuiAT + (ft)' ''■^^*-*<'^> ■ P**) 



IV. CONTACT TERMS AND ASYMPTOTICS OF THE SPECTRAL DENSITY 

The spectral sum rule ^ contains the quantity linii^^oo AG{uj, T). In this section we study the 
UV contact terms that appear in AG(a;, T), since those are the only ones that survive in the infinite 
frequency limit. We are working in infinite volume and start by considering the contribution of the 
energy density correlator. At finite separation in time, this correlator has a simple time-independent 
expression, 

j d^x (Too(x) Too(O)) = T^c, Vxq / . (39) 

The lattice sum rules, on the other hand, determine the correlator of Too integrated over all times. 
Combining these two pieces of information, one finds the contact term (see 16|, section 3.3), 

hm / d^xe''^^°(Too(x) Too(0))t-o = (-1^ + finite") (e - 3p) - -(e+p) , (40) 

where equation P^ has been used. 

We now coine to the contribution of the shear-stress correlator to AG. At large frequencies. 



the operator-product expansion (OPE) is applicable. The leading-order result is [19|] 26 1 

j d4xe--o(ri2(x)ri2(o))T-o "~°° -^(^ - 3^) - ^(^ + ^) • (41) 

Thus we see that 

AG(a;, T) "^-"^ 0(1) • (e - 3p) . (42) 

In particular, AG(a;,T) is finite in the w — >■ oo limit, and its value comes entirely from the 
contribution of the T^^j operator. 

The fact that both the zero-frequency correlator AG(0,T) and the contact term 
linitj^oo AG(a;, T) are finite implies that the correlator in coordinate space (i.e. as a function 
of Euclidean time t) diverges at short distances at most as 

AG(t, T) ~ ^^iilM . (linear combination of e, p) (43) 

at small t. From the spectral representation of this correlator 

/>oo 

AG(t, T) *^° / do; (Api2,i2(a;, T) - |Apoo,oo(^, T)) e""* , (44) 

it follows that the leading behavior of the shear spectral function is at most 

Api2,i2{'^,T) ^ rJ^ o-si^) ■ (linear combination of e, p) . (45) 



The leading order OPE result (j4ip does not allow us to obtain the finite part of the Wilson 
coefficient of the operator T^^ in the OPE of linii^^Qo AG {uj,T). However, the lattice sum rules 



3tain the finite part of the {T^^)t-o = e — 3p coefficient in the OPE of the energy 



allow us to o 
correlator J27l ]. 

hm / d4xe^'^^o(roo(x)Too(0))T-o = ce(e - 3p) - (f + 0{gl)) (e+p). (46) 

Indeed, the coefficient takes the value [la] 

^0 d log a 

where Eq. (|19p has been used in the second equality. As shown above, the leading term, which 
diverges logarithmically in the lattice spacing, cancels in AG(a;,T). The 0(5'g) terms vanish in the 
continuum limit. For Nc = 3, the finite terms in Eq. (|47p numerically amount to, respectively, 

lim {ce + - — ^ I = 2.3942 + 0.4215 - 0.2500 . (48) 

90^0 V 2boglJ 

The dominant contribution comes from the regularization-dependent {d^ — c'^) term. It is associated 

with the lack of continuous translation invariance on the lattice. We expect the regularization 

dependence to cancel in AG(cj,r), since it is a finite correlator. Therefore we take the two other 

terms to be representative of the size of T^^'s Wilson coefficient in the OPE of AG{uj, T) (after the 

inclusion of a factor |, see Eq. ([3T]) ). In the SU(3) gau ge t heory for T > 2.5Tc, (e — 3p)/T^ < 0.5 



while e/r^ w 4.5 is very weakly temperature-dependent 20]. From the size of the Wilson coefficient 
cg, we therefore expect Wm^^oo AG{uj,T) in the shear sum rule Eq. ([T]) to represent at most a 5% 
correction to the right-hand side of the equation. 

V. FINAL COMMENTS 

It is hoped that the spectral sum rule ([I]) will be useful in constraining the spectral density 
at finite temperature, especially when combined with lattice Monte-Carlo data on the Euclidean 
correlator. In order to be operational at all temperatures, the remaining contribution to AG{uj,T) 
should be computed in the large-frequency limit. Since we have seen that one contribution to AG 
has a dependence on the regularization, it is important to compute AG as a whole in the same 
regularization. Dimensional regularization is then probably the computationally most economic 
choice. 



The sum rule ([I]) is particularly useful in t 



tensor channel is very large, m2/Tc w 7.9 



re pure 



21 



gauge theory because the mass gap in the 



2^[28|. Thus even at T = 3Tc, the substracted 



10 

spectral density appearing in ([1]) coincides with the spectral density up to a; ~ 2.6T, in particular 
the latter is unaffected by the subtraction in the range of frequencies that determine the transport 
properties. Moreover, the contribution of the lightest tensor glueball to /9i2,i2(w, 0, T = 0) has 



u 



23( 1 ■ It is parametrized by a matrix element Ft in the notation 



been calculated on the lattice 

of 23(. At 3Tc, this contribution to the spectral sum rule ([TJ is small, 



Pi2,i2{u:, 0, 0) = ,7,^,, « 2%. (49) 



M3Z^W^,-e c.^ ■ - - m2 M|M 



At high temperatures, there is another contribution to the spectral sum rule that can be clearly 
isolated and identified. For a weakly coupled field theory, where the spectral density admits a 
transport peak at the origin of width much smaller than temperature, Teaney derived a sum 



rule [2J], 



r-A 



/A J 1 

— pi2,i2{u:,p,T)^-{e+p){vl). (50) 

-A ^ O 

Here A serves as a separation between the transport scale and the thermal scale. For ultra- 

relativistic quasiparticles, as one might expect to find in the high-temperature gluon plasma, {v^) = 

1. In that case, comparison of (j50p with ([1]) reveals that at high temperatures, the spectral integral 

([1]) receives an ~ 40% contribution from the transport peak. This means in particular that the 

prospect of determining the area under the transport peak from a numerically determined Euclidean 

correlator is realistic. 
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